The Predicted Signature of Neutrino Emission in Observations of 

Pulsating Pre- White Dwarf Stars 
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ABSTRACT 

Pre-white dwarf (PWD) evolution can be driven by energy losses from neutrino interactions 
in the core. Unlike solar neutrinos, these are not the by-product of nuclear fusion, but instead 
result from electron scattering processes in the hot, dense regions of the PWD core. 

We show that the observed rate of period change in cool PWD pulsators will constrain neu- 
trino emission in their cores, and we identify appropriate targets for future observation. Such a 
measurement will tell us whether the theories of lepton interactions correctly describe the pro- 
duction rates and therefore neutrino cooling of PWD evolution. This would represent the first 
test of standard lepton theory in dense plasma. 



Subject headings: dense matter — elementary particles — plasmas — stars: interiors 
GW Virginis — white dwarfs 



stars: variables: 



1. Introduction 

In general, stars are too remote — and observ- 
ables too few — to make them practical experimen- 
tal physics test-beds: our data are spent in sim- 
ply describing the dimensions of the objects under 
study. In many cases we must extrapolate experi- 
mental data over many orders of magnitude, or re- 
sort to untested calculations from first principles, 
to reach the regions of phase space that apply to 
stellar interiors. 

If we hope to overcome these problems and pur- 
sue "experimental" astrophysics, we can either at- 
tempt to increase the number of observables or 
find simpler stars. As first realized by Mestel 
(1952), the evolution of white dwarfs and pre- 
white dwarfs (PWDs) is primarily a simple cooling 
problem. In general, our growing understanding of 
white dwarf interiors and evolution has paralleled 
advances in the theory of dense plasmas, with the 
recognition of important influences like electron 
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degeneracy (Chandrasekar 1939), Coulomb inter- 
actions (Salpeter 1961), crystallization (Kirzhnitz 
1960; Abrikosov 1960; Salpeter 1961; Stevenson 
1980) and neutrino cooling effects (Chin, Chiu, 
& Stothers 1966; Winget, Hansen, & Van Horn 
1983; Kawaler, Hansen, & Winget 1985). Iben 
& Tutukov (1984) summarize the various mecha- 
nisms which dominate white dwarf evolution from 
the planetary nebula nucleus (PNN) stage to the 
coolest white dwarfs. 

On the observational side, the discovery of 
white dwarf pulsation in the 1960s, and pre-white 
dwarf pulsation in the 1970s, greatly increased 
the observable parameters available for compar- 
ison with theoretical models. These are short 
period, multiperiodic, g-mode variables, showing 
anywhere from a few to over a hundred separate 
periodicities on timescales of 100-3000 s. The pul- 
sating PWD stars are divided into two classes: the 
planetary nebula nucleus variables (PNNV stars), 
and the slightly more evolved GW Virginis (or 
simply GW Vir stars) which lack observed neb- 
ulae. With high surface gravities (log g ~ 6 — 7.5), 
and effective temperatures between 80,000 K and 
170,000 K, they occupy a region of the H-R dia- 
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gram between the high-T e ff end of the PNN branch 
and the top of the white dwarf cooling track. 
There are eight known PNNV stars, and four 
GW Vir stars (Ciardullo & Bond 1996). 

The evolutionary timescale of PWD stars is of 
order 10 6 years. During this short transition from 
PNN star to hot white dwarf, stellar radius and 
photon luminosity decrease by one and three or- 
ders of magnitude, respectively. High core den- 
sity and temperature allow electron scattering pro- 
cesses to produce a large neutrino flux which re- 
mains roughly constant during this time. As pho- 
ton luminosity plummets, neutrinos contribute an 
increasing fraction of the total energy losses. Neu- 
trino emission eventually comes to dominate the 
overall evolution of the star. 

Unlike photon energy which must diffuse rela- 
tively slowly through the entire star before emerg- 
ing into space, neutrinos created near the center 
of the PWD escape directly. This neutrino lumi- 
nosity cools the center of the star, maintaining a 
temperature inversion similar to that within stars 
at the tip of the red giant branch. Calculations of 
the relevant reaction rates were performed initially 
by Beaudet, Petrosian, & Salpeter (1967) based on 
the theory of weak interactions proposed by Feyn- 
man & Gell-Mann (1958). Later, Dicus (1972) 
and Dicus et al. (1976) recalculated these rates 
in the unified electroweak theory of Weinberg and 
Salam (Weinberg 1967, Salam 1968). All of these 
calculations are theoretical, however. We have no 
direct experimental or observational confirmation 
of neutrino production rates under conditions ap- 
propriate to PWD interiors. 

The cooling of a GW Vir interior tends to in- 
crease the periods of each given pulsation mode. 
Their high luminosity (log L ~ — 3) means they 
cool much more rapidly than cooler white dwarf 
variables. GW Vir period changes are therefore 
expected to be more rapid also. Winget, Hansen 
& Van Horn (1983) show that the e-folding time 
for period changes in GW Vir stars should be of 
the same order as the evolutionary timescale — 
10 6 years; such rapid changes are measurable in 
1 — 3 years time. This is an exciting prospect: to 
measure directly, on human timescales, the rate 
of evolution of a star, and specifically to place 
strict constraints on the mechanisms which reg- 
ulate the evolution of a stellar interior. Over 30 
years ago, Chin, Chiu, and Stothers (1966) pre- 



dicted that at some point in PWD evolution neu- 
trino losses should dominate all other cooling pro- 
cesses. Asteroseismological analysis can tell us 
which stars these are, and then measurement of 
period changes can tell us if our neutrino physics 
is right. 

Such a test has implications far beyond the 
study of PWD evolution. For instance, one of 
the fundamental questions of stellar astrophysics 
is the length of time stars spend on the main 
sequence. Answering this question requires pre- 
cise knowledge of the p-p and CNO nuclear re- 
action rates. Currently, the best laboratory for 
measuring these rates is our own Sun, since ter- 
restrial labs cannot in general reproduce the con- 
ditions of the stellar interior. However, models 
which successfully reproduce the known structure 
of the Sun predict a neutrino flux two to three 
times that measured by earthly detectors (Bah- 
call & Pinsonneault 1996, and references therein). 
For a long time, it was thought the problem might 
reside in our incomplete knowledge of conditions 
in the solar interior. Recently, helioseismology 
projects such as the Global Oscillation Network 
Group (GONG) have resulted in the measurement 
of millions of solar pulsation frequencies (Harvey 
et al. 1996). With so many parameters to con- 
strain model properties, the possibility that the 
solar neutrino problem can be solved through vari- 
ations in the thermodynamics or mechanics seems 
to be excluded (Bahcall & Pinsonneault 1996). 
The problem, then, almost certainly lies with the 
way we handle the nuclear physics. 

Under the most intense scrutiny is the stan- 
dard theory of lepton interactions. Our calcula- 
tions of neutrino emission from PWDs are based 
on this same theory. In PWDs, however, the en- 
ergy loss rate due to neutrinos is thousands of 
times greater than in the Sun. Measurement of 
the effects of neutrino interactions in PWDs would 
afford a critical independent test not only of the 
standard lepton theory but also of non-standard 
theories brought forward to solve the solar neu- 
trino problem. 

To explore this possibility, we calculated PWD 
evolutionary tracks using different neutrino pro- 
duction rates. In the next section we describe 
the calculation of those rates and summarize the 
basic interactions that lead to neutrino emission 
in PWD interiors. Section 3 describes PWD se- 
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qucnces with varied neutrino production rates and 
examines effects on measurable quantities such as 
T e ff, surface gravity, and rate of period change. 
Finally, in § 4 we discuss prospects for placing ob- 
servational constraints on neutrino physics, and 
we identify appropriate targets for future observa- 
tion. 

2. Neutrino Cooling in Pre- White Dwarf 
Interiors 

Unlike the solar neutrino flux, neutrino emis- 
sion in PWDs is not a by-product of nuclear fu- 
sion. Instead, the density and temperature in their 
cores are high enough (log p c ~ 6-7, log T c ~ 7-8) 
to produce neutrinos directly through several dif- 
ferent scattering processes. The two most impor- 
tant processes are neutrino bremmstrahlung and 
plasmon excitation. Neutrino bremmstrahlung 
is much like the normal bremmstrahlung process 
familiar to astrophysicists, in which high-energy 
electrons scatter off nuclei, emitting X-rays. At 
the high density and temperature of PWD in- 
teriors, however, neutrinos can be produced in- 
stead. These same conditions support the ex- 
istence of thermally excited photons within the 
plasma, analogous to phonons propagating within 
a metal lattice. These "plasmons" have a finite 
lifetime and decay to form a neutrino and antineu- 
trino. 3 

The possible relevance of the plasma process 
to stellar astrophysics was first pointed out by 
Adams, Rudcrman, & Woo (1963), who subse- 
quently calculated rates based on the theory of 
Feynman & Gell-Mann (1958). Beaudet, Pet- 
rosian & Salpeter (1967) were the first to incor- 
porate them into stellar evolution calculations. 
Later, Dicus (1972) recalculated the rates of vari- 
ous neutrino processes in the unified electro-weak 
theory of Weinberg and Salam (Weinberg 1967, 
Salam 1968). 

The rates used in our stellar evolution code, 
ISUEVO, derive from updated calculations by Itoh 
et al. (1996), and include the plasmon, bremm- 

3 Actually there are two types of plasmons. The process 
described here is that of the transverse plasmon. The other, 
longitudinal plasmon corresponds to an oscillation in the 
electron gas similar to a sound wave, but is usually less 
important as a neutrino source in hot white dwarfs (Itoh 
ct al. 1992). 



strahlung, and several less important neutrino pro- 
duction processes. The evolution code ISUEVO 
(Dehner 1996; see also Dehner & Kawaler 1995) is 
optimized for the construction of PWD and white 
dwarf models. The models used in this investiga- 
tion are based on the evolution of a 3 M Q model 
from the Zero Age Main Sequence through the 
thermally pulsing AGB phase. After reaching a 
stable thermally pulsing stage (about 15 thermal 
pulses), mass loss was invoked until the model 
evolved to high temperatures. This model (rep- 
resenting a PNN) had a final mass of ~ 0.6 M©, 
and a helium-rich outer layer. Additional details 
concerning the construction of this evolution se- 
quence (and others of different mass, discussed in 
§ 3, below) can be found in O'Brien (1998). 

To study the direct effects of neutrino losses on 
PWD evolution, we introduced artificially altered 
rates well before the evolving models reached the 
PWD track. If we simply changed the rates be- 
ginning at the hot end of the PWD sequence, the 
thermal structure of each model would take several 
thermal timescales to relax to a new equilibrium 
configuration based on the new rates. Unfortu- 
nately, this relaxation time is of the same order as 
the PWD cooling time, and so only the cool end 
of the sequence would see the full effects of the 
new rates on their evolutionary timescales. There- 
fore, the enhanced and diminished rates described 
in the next section were introduced into evolution- 
ary calculations beginning at the base of the AGB. 
The resulting thermal structure of the initial PWD 
"seed" models was then already consistent with 
the neutrino rates used during the prior evolution 
that produced them. 

3. Pre- White Dwarf Sequences with Dif- 
ferent Neutrino Rates 

Starting with the PWD seed models above, we 
evolved the models from high L and T e ff down 
toward the white dwarf cooling track. Three se- 
quences were calculated. The base sequence used 
the normal neutrino production rates. Another se- 
quence used rates diminished by a factor of three 
(at any given p and T in the stellar interior) over 
the normal rates, while the third sequence used 
rates enhanced by a factor of three. This trio spans 
nearly one order of magnitude in neutrino produc- 
tion. 
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The resulting 0.6 M Q evolutionary sequences 
are shown in Figure 1, from T e g ~ 170, 000 if— 
equivalent to the hottest PWDs known — down to 
about 35, 000 if. Luminosity decreases by almost 
four orders of magnitude in approximately five 
million years. The GW Vir instability strip oc- 
cupies the left half of the figure, above T e g ~ 
80, 000 if (log T e g = 4.9), a temperature reached 
by the PWD models in only 500,000 years. 
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Fig. 1. — Evolutionary tracks for three 0.6 M 
model sequences with different input neutrino pro- 
duction rates. The upper and lower tracks were 
calculated with rates one-third and three times the 
normal rates (middle track), respectively. 

The most striking aspect of Figure 1 is the sim- 
ilarity of the tracks: changing the neutrino rates 
seems to have little effect on the luminosity at 
a given T c g at any point in PWD evolution, de- 
spite the importance of neutrino losses as a cool- 
ing mechanism over much of this range. In Fig- 
ure 2, we find that, for all three sequences, neu- 
trino losses are the primary cooling mechanism 
over the approximate range 100, 000 if < T c g < 
30, 000 if. Plasmon reactions dominate over the 
bremmstrahlung process for 0.6 M Q models at all 
stages of PWD evolution, as shown in Figure 3. 

The ratio L v jL 1 also increases with stellar 
mass. In the T cS range 80,000-100,000 K, L v /L 1 
for a 0.66 M Q model sequence is nearly 30% higher 
than for a 0.60 M Q sequence. 

Figures 1 and 2 show that the differences in L 
and T e ff are smallest when the neutrinos are im- 
portant. This is because the primary structural 
effect of changing the neutrino rates is on the ra- 



dius of the models (Figure 4), causing the tracks 
to assume a position in the L-T e g plane normally 
occupied by models of slightly higher mass (for en- 
hanced rates) or lower mass (for diminished rates) . 
However, at lower temperatures electron degener- 
acy becomes increasingly important as a mechan- 
ical support against gravity (and thus in deter- 
mining the final stellar radius); neutrino cooling 
only affects the thermal processes participating in 
the mechanical structure. Even at high luminos- 
ity, however, different neutrino rates result in only 
small changes in measurable quantities such as 
surface gravity. Current observational techniques 
could not hope to resolve such small differences. 
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Fig. 2. — Ratio of the neutrino luminosity to 
the photon luminosity as a function of T e g, for 
three 0.6 Mq model sequences with different input 
neutrino production rates. 

Figure 5 shows a more tangible effect of chang- 
ing the rates. Even though models with differ- 
ent rates look much the same at a given T c ff, 
they get there at widely differing times, since the 
rate of evolution along a track is directly depen- 
dent on the importance of neutrino emission as a 
source of cooling. For example, the model with en- 
hanced neutrino rates cools from 100, 000 K down 
to 65,000 K in 600,000 years, while the model 
with diminished neutrino rates takes 1.3 million 
years, more than twice as long, to cool by the 
same amount. The maximum difference in the 
slope of the different curves in Figure 5 occurs at 
T cS ~ 80, 000 K. Thus the epoch where the rate 
of evolution is most sensitive to the assumed neu- 
trino physics corresponds to the position in the 
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H-R diagram occupied by the coolest pulsators in 
the PWD instability strip. On the other hand, for 
stars in the strip hotter than 100, 000 K Figure 5 
shows that evolutionary rates do not depend on 
neutrino rates. 
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to that of the next. To determine dll/dt, a bet- 
ter technique is to measure the cumulative phase 
change in a mode with a slowly changing period. 
This is accomplished by comparing the observed 
times of maxima (O) in the light curve to the times 
of maxima (C) calculated from an assumption of 
constant period. The resulting plot of (O — C) 
shows the phase drift associated with a changing 
period. A constant rate of period change, dll/dt, 
enters as a quadratic term in time: 

where U to is the period at time t (see for example 
Winget et al. 1985, 1991 and Kepler et al. 1995). 
To measure dll/dt with confidence, the star must 
of course have stable and fully resolved pulsation 
periods, with reliable phase measurements from 
season to season. 



Fig. 3. — Fraction of total neutrino lu- 
minosity contributed by the plasmon and 
bremsstrahlung processes, as a function of T e g, for 
three 0.6 M© model sequences with different neu- 
trino production rates. 

Our expectations are borne out in Figure 6, 
which shows the rate of change in period, Yl/Yl 
(= d(lnU)/dt), as a function of period, II, for 
PWD models at 140,000 K (lower panel) and 
80,000 K (upper panel), given normal, enhanced, 
and diminished neutrino production rates. The 
rate of period change U/H in the cooler models 
changes by a factor of four between the enhanced 
and diminished rates. Changing the neutrino rates 
has little effect on II/II in the hotter model, con- 
sistent with the results from Figure 5. We now 
turn to the exciting implications of these results, 
and explore the possibility and practicality of mea- 
suring II/n in cool pulsating PWD stars. We can 
then identify likely targets for future observational 
campaigns. 

4. Prospects for Measuring Neutrino Cool- 
ing Effects 

4.1. Determination of dll/dt 

Unfortunately, the period changes expected to 
occur in PWD stars are too small to detect from 
simple comparison of the period from one year 
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Fig. 4. — Thermal structure at three different evo- 
lutionary stages, T e g — 140,000 K (upper panel), 
80,000 K (middle panel), and 30,000 K (lower 
panel), for 0.6 M models with different neutrino 
production rates. 
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Fig. 5. — Evolution of T e fj with time for 
three O.6M model sequences with different input 
neutrino production rates. 



4.2. Prospective Targets 

Measurements of secular period change, n/n, 
in white dwarfs have been attempted by a num- 
ber of investigations, with either measurements 
made or tight upper limits set for the GW Vir 
star PG 1159 (Winget et al. 1985, 1991, Costa & 
Kepler 1998) and G117-B15A (Kepler et al. 1995). 
Unfortunately, neutrino cooling is not expected to 
be an important effect for either of these stars. On 
the other hand PG 0122, with a T cff of 80,000 K, 
occupies the stage in GW Vir evolution most 
highly dominated by neutrino emission. O'Brien 
et al. (1998) show that PG 0122 is in addition the 
most massive GW Vir star, which should enhance 
neutrino effects as well. 



Kawaler, Hansen, & Winget (1985) and Kawaler 
& Bradley (1994) present predicted values of 
dll/dt for models relevant to GW Vir and PNNV 
stars; the only observed value of dn/dt, that for 
PG 1159 itself, is consistent with these models. 
However, as Kawaler & Bradley (1994) demon- 
strated, for a star as hot as PG 1159, dn/dt is 
strongly affected by mode trapping. This is an 
effect whereby some modes become excluded from 
regions below subsurface composition discontinu- 
ities. Kawaler & Bradley (1994) show that, in 
general, dn/dt should be positive; this reflects 
the overall cooling of the model (Winget, Hansen, 
& Van Horn 1983). Trapped modes, however, 
are concentrated in the outer layers, within which 
contraction dominates cooling; therefore trapped 
modes can show periods which decrease with time. 
Thus, mode trapping can complicate the interpre- 
tation of measured period changes in hot PWDs. 
As GW Vir stars cool, the surface contraction rate 
decreases relative to the cooling rate of the inte- 
rior. So, while mode trapping can still influence 
the pulsation period distribution itself, the rates of 
period change become more similar from mode to 
mode in cooler GW Vir stars. Kawaler & Bradley 
(1994) found that the sign of dn/dt could be dif- 
ferent for different modes in hot GW Vir models; 
by the time those models evolve to the cool end of 
the strip, the period change rates are all positive. 
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Fig. 6. — Rate of period change for 0.6 M models 
with different neutrino rates. 

In order to measure IT/n with confidence, a star 
must have a stable and fully resolved pulsation, 
with stable phase measurements from season to 
season. PG 0122 is a very stable pulsator: over the 
past decade, it has shown a consistent pulsation 
spectrum, with the large-amplitude modes present 
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at the same frequencies during each of three in- 
tensive observing seasons in 1986, 1990, and 1996. 
The amplitudes of each of the dominant modes re- 
mained approximately constant as well (O'Brien 
ct al. 1998). Therefore, PG 0122 is an excellent 
candidate for measurement of the rate of secular 
period change caused by the evolutionary cooling 
of its interior. 

In addition to the physics governing neutrino 
production, PG 0122 is an ideal target for measur- 
ing neutrino emission rates because of the minimal 
influence of any mode trapping on interpretation 
of its n/II. As mentioned above, for stars below 
100,000 K trapping no longer significantly affects 
IT/n from mode to mode. 

From Figure 6, we estimate the value of n/II 
for the dominant pulsation mode (II = 400 s) in 
PG 0122 to be about 6 x 10~ 15 sec - 1 . With 
this rate of period change, the period should in- 
crease by about 0.001 s in 10 years; this is smaller 
than the period uncertainty for a run length of 
several months (assuming a frequency precision of 
10 x run length )- However, the accumulated phase 
advance over a ten year period should be nearly 
two full cycles. 

Using the periods alone from the 1986, 1990, 
and 1996 data, O'Brien et al. (1998) attempted 
to calculate U/U directly. From the best least- 
squares periods from 1996 and 1986, they calcu- 
late a period change of —0.10 ± 0.02 s, implying 
ri/Tl = -7 x 10~ 13 sec _1 , about 100 times larger 
in magnitude — and different in sign — than theory 
expects. However, as O'Brien et al. (1998) point 
out, this calculation is based on the formal errors 
from a least-squares fit to the observed periods, 
and the formal least-squares error generally un- 
derestimates the true error by an order of mag- 
nitude. In practice, the data currently available 
allow an upper limit to be set on the absolute mag- 
nitude of Il/n for PG 0122 of 1.5 x 10~ 12 sec _1 . 
In view of the importance of measuring IT/TI for 
this star — as well as for the other cool GW Vir 
stars — we must continue analysis of archival data 
and mount observing campaigns in the near fu- 
ture to monitor all the known GW Vir stars with 
T c ff < 100, 000 K. With frequent observation, an 
accumulated phase advance of half a cycle, com- 
bined with the techniques described above, could 
be used to determine n/II for the GW Vir stars 
PG 0122, PG 2131, and PG 1707 in two to three 



years. In the case of PG 0122, the data presented 
in O'Brien ct al. (1998) provide a key anchor for 
this investigation. 

5. Summary and Conclusions 

We have shown that the predicted rates of pe- 
riod change in GW Vir stars near the cool end of 
the instability strip are sensitive to the neutrino 
production rates used in stellar models. The per- 
sistence of the solar neutrino problem has made 
the standard model of neutrino interactions one 
of the most intensely scrutinized theories in all of 
physics. Determination of II in the GW Vir stars 
PG 0122 and PG 2131 will provide an important 
test of the standard model and of any new theories 
put forward to replace it. 
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